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A method for incorporating systematic errors into branching ratio limits which are not obtained 
from a simple counting analysis has been suggested by Mark Convery [1]. The derivation makes 
some approximations which are not necessarily valid. This note presents the full solution as an 
alternative. The basic idea is a simple extension of the Cousins and Highland philosophy [2]. 
Before systematics are considered, an analysis using a maximum likelihood fit returns a central 
value for the branching ratio (B) and a statistical error The likelihood function is 
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Following the Convery notation, we associate 5* with the nominal efficiency and as as the (Gaus- 
sian) error on the efficiency. Adding the uncertainty on the efficiency changes the likelihood to: 
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From Mathematica®, the integral in Eq. 2 is: 
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Removing unimportant multiplicative constants and changing variables from as to a^ = as/ S 
gives: 
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It turns out that as long as the efficiency S is sufficiently small (generally less than 10% but 
dependent on other parameters), the second erf term evaluates to —1 and the dependence on the 
efficiency is removed. 

The solution to the integral presented by Convery (for as ^ S) can be written as: 
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The differences between Eq. 3 and Eq. 4 are the two erf terms in Eq. 3. The first erf term 
affects the tails of the distribution and becomes increasingly important as a^ increases. The second 
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erf term affects the peak position and is important when S" ± crs' is not easily contained in the 
region {0, 1}. Or, for a fixed 0"^, when S approaches unity. Next we compare the two results after 
modifying Equations 4 and 3 to normalize them such that p{B = B) = 1. 

First we check the effect for relatively large and small S for which the first erf term becomes 
important. Each plot of Figure 1 shows a comparison between the full solution in red and the 
approximate solution in black. There is very little discernible difference between the two solutions. 
The different plots show results for B = 0.5, B = —0.5, and B = —1.5. To set an upper limit, one 
often integrates the probability over the physical region only {B > 0). Figure 2 shows the results 
for p{B) over the range B G {0, 17} for the case of 5 = —1.5 and as = 0.5 which corresponds 
to a 3a negative fluctuation. In this case clear differences between the full solution (in red) and 
the approximate solution (black) can be seen for > 0.3. Note that Fig. 1(a) and Fig. 2(b) 
show the same curves, only the range has changed. Clearly an attempt to flnd an upper limit 
by integrating the area under the approximate solution is problematic for all the cases shown in 
Fig. 2. Conversely, the full solution flnds an acceptable upper limit. 
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Figure 1: Each plot shows a comparison of the approximate solution given by Eq. 4 in black to 
the full solution given by Eq. 3 in red. For all plots, cb = 0.5, = 0.4, and S = e = 0.1. The 
three plots show results for B = 0.5, B = —0.5, and B = —1.5. 



Second we check the effect of the second erf term of Eq. 3 which is important when the 
integration of efficiency from to 1 in Eq. 2 cuts off a significant part of the Gaussian defined by 
S ± as = S ± a^S. Figure 3(a) is a repeat of Fig. 1(a) on a different scale and again shows little 
difference between the two methods. Figures 3(b) and 3(c) show the effect of the second erf term 
as S* ^ 1. 

In conclusion, Eq. 3 provides a more exact and robust implementation of the original suggestion 
by Convery [1] on incorporating multiplicative systematic uncertainties in branching ratio limits. 
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Figure 2: Each plot shows a comparison of the approximate solution given by Eq. 4 in black to 
the full solution given by Eq. 3 in red. For all plots, 0"^ = 0.5, B = —1.5, and 5" = e = 0.1. The 
three plots show results for = 0.3, = 0.4, and 0"^ = 0.5. In this case, the full solution is 
indistinguishable from the full solution with the second erf term replaced by —1. 




Figure 3: Each plot shows a comparison of the approximate solution given by Eq. 4 in black to 
the full solution given by Eq. 3 in red. For all plots, = 0.5, B = 0.5, and = 0.4. The three 
plots show results for S" = e = 0.1, S = e = 0.8, and S = e = 0.9. In this case, the full solution is 
nearly indistinguishable from the full solution with the first erf term replaced by +1. 
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